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We discuss the effect of kinetic energy of the relative motion becoming spurious for 
separate fragments on the selfconsistent mean-field fission barriers. The treatment of 
the relative motion in the cluster model is contrasted with the necessity of a simpler and 
approximate approach in the mean-field theory. A scheme of the energy correction to 
the Hartree-Fock is proposed. The results obtained with the effective Skyrme interaction 
SLy6 show that the correction, previously estimated as ~ 8 MeV in A = 70 — 100 nuclei, 
amounts to 4 MeV in the medium heavy nucleus 198 Hg and to null in 238 U. However, 
the corrected barrier implies a shorter fission half-life of the latter nucleus. The same 
effect is expected to lower barriers for multipartition (i.e. ternary fission, etc) and make 
hyperdeformed minima less stable. 



1. Introduction 

It seems that the existing calculations of fission barriers overestimate energies of 
nuclear configurations close to scission by including a spurious contribution of ki- 
netic energy of the fragments' relative motion. This may be seen as follows: The 
binding of two separated selfbound systems is equal to E sep = B\ + B2 + Vint, 
with Bi fragment binding energies and V, n t Coulomb energy. For a compound sys- 
tem, energy E(l + 2), should tend to E sep for separate entities, hence it should 
not contain fragment center of mass (cm.) energies E c . m \i). However, in a stan- 
dard Hartree-Fock (HF), energy of a two-piece configuration still contains the term 
E CAn .(rel) — E r r r , (1) +E r .. m _ (2) — E c , m _ (1 + 2) , corresponding to the relative motion 
of two fragments 

Within the mean-field theory, this overestimate arises from the cm. kinetic en- 
ergy correction that has to be subtracted from the expectation value of the Hamil- 
tonian (energy functional) to obtain the binding. The expectation value of the 
operator (J2 k p k ) 2 /(2AM) on the Slater state reads: 

E c . m .{A) = -i- |X> I p 2 I k) -J2 I (* I P 1 1 2 1 . (1) 



,fc=l k^l 



1 
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with k, I labelling occupied single particle states, and M nucleon mass. One should 
include one ~E c . m . correction for a compound system, two corrections for two sepa- 
rate fragments, three corrections for three fragments, etc. For two distant fragments, 
A1+A2 = A, if one could distinguish particles belonging to each, which would imply 
the vanishing expectation value ((J^feei P fe ) ' Cl2i£2Pl))i 

E c . m .{A x ) + E c . m .(A 2 ) = E c . m .(A) + ^^M+AWM , (2) 

The second term on the right-hand side is just the asymptotic (for the distance 
R — ► 00) value of kinetic energy of the fragments' relative motion, E c _ m .(rel) =< 
P 2 rel > /(2/x),with P rPj = (A 2 /A)Y: kel Pk-{Ai/A)J2ke2Pk, and n = MA-lA 2 /A. 
This quantity becomes spurious for a two-piece system as it does not contribute 
to the binding. As noted by Berger and Gogny, ^ this asymptotic term should be 
subtracted from the HF energy to obtain proper fusion barriers. 

In practical HF calculations the -E c . m . correction is included in various forms 
(see^. Here wc confine the discussion to effective interactions which use the natural 
definition (JTJ) . With such interactions one obtains E c , m , — 5.5 — 8 MeV for A = 
40 — 250, decreasing with A. This quantity should be subtracted from E(l + 2) 
in a consistent theory: partially - for shapes with constriction, totally - for two- 
piece systems. This subtraction is usually included in calculations of fusion barriers 
(otherwise the barriers are too highE| , ), but omitted in fission studies. However, 
even for fusion barriers, a gradual dependence of the correction on the compactness 
of the system is missing. 

As found in Ref. the subtraction of the asymptotic value of E c , m (rel) brings 
the calculated HF fission barriers closer to the experimental values in medium-size 
A = 70 — 100 nuclei. In the present work we estimate the effect of the shape- 
dependent correction E c/m ,{rel,shape)on fission barriers in heavier nuclei, using 
the Skyrme effective interaction SLy6^ (section 3). The correction is discussed and 
defined in section 2, where we also consider a different treatment of the relative 
kinetic energy in the cluster model and in the mean-field theory. Conclusions are 
given in section 4. 

It is remarkable that a correction of similar property and magnitude, although 
based on completely different grounds, has been introduced in macroscopic- 
microscopic 

calculations™ in rder to obtain a better agreement with the exper- 
imenal fission barriers in the same regions of nuclei. 

2. General discussion 

The difficulty in the determination of the fission or fusion barrier discussed here is 
pertinent to nuclear models in which the binding energy of the far separated nuclei 
1 and 2 treated as a one system is different from B\ + B 2 , with Bi determined 
separately. Since the standard HF belongs to this category, it needs a correction 
which would ensure that £"(1 + 2) — > B\ + B 2 in fission (with E(l + 2) understood 
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as adiabatic energy) . The smoothness of such correction with the evolving nuclear 
shape is a natural requirement. 

At the heart of the difficulty lies the identity of particles that impedes a 
definition of the relative coordinate and motion of two interacting subsystems. 
This problem is crucial in studies of light nuclei, where it is treated within 
the resonatig group method (RGM) that is basically a version of the genera- 
tor coordinate method (GCM). A cluster configuration A\ + A2 is represented 
as an expansion $a = J dr<p(r)$ T onto an overcomplete basis formed by states 

= AaiA 2 [<K r ~ *rel)®Ai$A 3 ] with $a« the cluster states depending on 

intrinsic coordinates, ip the amplitude of the relative motion and AaxA 2 the anti- 
symmctrizer containing permutations mixing the coordinates of the first A\ with 
those of the last A2 nucleons. Thus, the label r of the basis states assumes the role 
of the intercluster coordinate. 

From the Schrodinger equation for $a one obtains the Hill- Wheeler equation 
for the amplitude of the relative motion in coordinate r, with a well defined Hamil- 
tonian. However, a decomposition of this Hamiltonian into kinetic and potential 
parts, T + U, is arbitrary: The relative kinetic energy operator is assu med bs 
T = — (fi,V r ) 2 /(2/i) with the reduced mass fi, and this fixes the potential Z//Ep^l. 
So obtained potentials are much deeper in the compound nucleus region than those 
implied by the mean field; in this way the Pauli exclusion influences the relative 
motion of the overlaping clusters. 

It follows that while the RGM (or GCM) provides a solution to the rela- 
tive motion problem, its ingredients, like the potential U, do not seem to be of 
much use for the mean-field theory. As an aplication of the full RGM (GCM) 
method for heavy nuclei would be prohibitively difficult, one would rather include 
a kinetic energy correction in the relatively simple HF method to improve en- 
ergy asymptotics for separated clusters. However, this cannot be just the expecta- 
tion value of [A/ (2M AiA2)]P 2 rel in the Slater state: Owing to the incompatibility 
of the P r ei variable with the antisymmetry of the Slater determinant, its value, 
(< p 2 > + I ( fc I P I I 2 l\A{A - 1)])/(2M), is by ~ 10 MeV larger than the 

proper value of E c , m .(rel). This is why obtaining the correct value of Eq.([2]) requires 
an extention that goes beyond HF. Some guidance might be provided by realistic 
internuclear potentials used in fusion studies, e.g. ^1 These potentials have cor- 
rect two-fragments asymptotics, but should be smoothly replaced by the mean-field 
energy for compact mononuclear configurations, where they become irrelevant. 

One possibility to proceed^ is to introduce a measure of the fragment separation 
£ which would replace the relative distance r and define the subtracted portion of 
the relative kinetic energy, E c _ m _(r el, shape) = £E c , m .(rel). To this aim, consider 
dispersion of the number of particles in the fc-th HF orbital, residing in the volume 
Vi of the first fragment with A\ nucleons. This reads Pk(l—Pk), withpfc = Jy | ipk | 2 
(with many-particle correlations ignored). For completely divided fragments the k- 
th wave function is localized, so pk =0 or 1 and dispersion vanishes. We define £ by 
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means of dispersion averaged over the occupied orbitals 

, , |'2\ Ek-occPk(l ~ Pk) f 2 \ Ek-occ(Pk -P? () 

5 U; p(i-p) w p(i-p) ' lj 

with p = 2 X)fc-occP fc /^ = A±/A (with obvious modifications for pairing included). 
So defined £ is positive, reaches the maximal value 1 for separated fragments and 
falls to zero for wave functions uniformly smeared over two fragments. It has thus 
necessary properties to show main effects of a gradual inclusion of E c _ m .(rel) in HF 
energy. 

Subtraction of a varying portion of E c . m .(rel) will change energy balance be- 
tween configurations with and without constriction, lowering the former with re- 
spect to the latter. As found in the study ^ of fusion barriers, values of £ at the 
barrier vary between 0.7 and 0.9 and decrease together with the interfragment 
distance R. The latter is defined as the distance between the centers of mass of 
two half-spaces, containing A\ and A 2 nucleons. Clearly, fusion barriers calculated 
using the shape-dependent correction £E c . m _(rel) are higher than those obtained 
by subtracting the whole asymptotic value E c _ m _(rel). The related increase in the 
barrier height will depend on the slope dV/dR (without any correction): a small 
increase for a large positive slope, a larger increase (and the barrier shift towards 
smaller R) for small positive or negative slopes. For example, with the SLy6 force, 
the inclusion of the ^-dependence leads to the increase in fusion barrier by 1.8 MeV 
for 48 Ca+ 48 Ca and by 2.5 MeV MeV for 48 Ca+ 208 Pb, with the inward shift of the 
barrier top by 0.6 and 1.5 fm, respectively^. 

In the present study of fission barriers we use a different prescription for the 
relative kinetic energy correction. For a system of A nucleons consider its division 
A\ + A 2 into volumes V\ and V%. Calculate the quantity pk for each wave function 
and call it localized in V\ (V2) if Pk > 1 — e (pk < e), with some small e (we use 
e = 0.03). Suppose that for a given nuclear shape (configuration) N\ wave functions 
are localized in V\ and N 2 in V 2 . Then the correction for this shape is defined as 

E c . m .(rel, shape) = -^E^Nx) + -±E c . m .{N 2 ) - £ c .m.(iVi + N 2 ). (4) 

A\ A 2 A 

This quantity tends to for no localization and to the asymptotic value E c , m ,(rel) 
for divided fragments (full localization). The correction ^ is more directly related 
to the localized orbitals than £E c . m .(rel), but still not completely satisfactory. Ul- 
timately, it would be desirable to define the correction for relative kinetic energy 
as a part of the energy functional and treat it variationally. 

3. Shape-dependent correction to barriers 

We have calculated the fission barriers with and without the E c , m .(rel, shape) cor- 
rection in 198 Hg and 238 U. Pairing was included as the delta interaction in the BCS 
scheme, using the cutof according to the prescription of Ref. The delta interac- 
tion strength was fixed at V n = 316 MeV fm 3 for neutrons and V p = 322 MeV fm 3 
for protons. 
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Fig. 1. Fission barrier in 198 Hg without (squares) and with (circles) the E c . m .(rel, shape) correc- 
tion. 



The calculated fission barrier in 198 Hg (Fig. 1) is mass-symmetric (A\ = A2) 
and the saddle corresponds to a large elongation with the quadrupole moment close 
to Q = 300 b, cf Fig. 2. To relate Fig. 1 to other studies, we mention that the energy 
plot for Q < 20 b, calculated with the Gogny interaction, may be found in Fig. 2 of 

while the whole macroscopic-microscopic fission barrier was given in Fig. 3 of 
L^2J Between the secondary minimum at Q = 45 b and Q ~ 100 b the HF minima are 
soft to mass-asymmetry or even slightly mass- asymmetric. Our calculated barrier of 
27.5 MeV is lowered owing to the E c m (rel, shape) correction to 23.7 MeV vs. the 
experimental value of 20.4 MeV and 19.3 MeV calculated in^l It may be seen 
that the relative kinetic energy correction of 3.8 MeV at the 198 Hg fission barrier 
is smaller than those for A = 70 — 100 nuclei^, but still significant. 

The calculated barrier in 238 TJ is shown in Fig. 3. The triaxial first hump of 7.55 
MeV at Q — 60 b and the mass-asymmetric second hump of 8 MeV at Q = 125 
b (that agrees with the result with SLy6 reported in are both larger than 
the experimental barriers (inner and outer) close to 6 MeV 120213 At the second 
barrier, the neck is still not yet developed (cf shapes in Fig. 4) and single-particle 
orbitals are not well localized in parts of the volume corresponding to A\ = 138 
and A2 = 100, the partition chosen close to the maximum of the experimental mass 
yield 122. Hence the E c<m ,(rel, shape) correction vanishes and does not influence the 
barrier height in this case. For quadrupole moments Q > 125 b the localization 
begins and the correction grows slowly with Q. At Q = 200 b (shape in Fig. 4) it 
amounts already to 1.9 MeV. This implies that the barrier relevant for the quantum 
tunneling, while not being higher, becomes shorter. 

An estimate of the effect on the fission half-life may be obtained by using the 
relation A log T s f w 0.8686A5, with action S = J y/2B ef f(E-E g . s )d(3 2 , see e.g. 

For B e ff we can use the cranking mass parameters, typical for the appropriate 
deformation range, calculated with the Woods-Saxon potential. An approximate 
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Fig. 2. Barrier configuration in 198 Hg, Q =300 b; density contours 0.01 fm 3 apart. 



relation between the quadrupole moment and the deformation parameter /3 2 , #2 ~ 
\/5^Q/(3r^ 5 / 3 ), with the assumed r Q = 1.2 fm, gives /3 2 « 0.01 x Q for A = 238. 
Between /3 2 = 1.8 and 2.3, where the correction causes an appreciable difference 
in the integrand of action, the mass parameter B22 decreases from about 20 to 5 
H /MeV and is similar to B23, while -B33 increases from about 30 to 60 fi. 2 /MeV 
(subscripts of the mass tensor refer to deformation parameters /3a, ^j). Taking 
B e ff = 10h 2 /MeV (as /3 3 increases with (3 2 , Beff > B22), and estimating AS as 
(2.3 - 1.8) x y/2B ef fE av with E av equal to 1.0 MeV (cf Fig.3), we obtain a rough 
estimate AS = 2.23 and AlogT s / = 1.94. Thus, the expected change in the fission 
half-life is abo ut tw o orders of magnitude, compared to the experimental value of 
logT s f of 23.4 IS. The true correction AlogT s / is probably at least that large: 
The recently calculated selfconsistent cranking inertia parameters 1221 that should 
be used in the exact calculation, seem larger than the Woods-Saxon cranking mass 
parameters. 

It is worth mentioning that the calculated barriers will be still lowered by the 
rotational correction. In 238 U, one can expect more than a 1 MeV correction at the 
first barrier , an d more than a 2 MeV correction at the second barrier, based on 
calculations ^1 An even larger rotational correction should be expected for 198 Hg 
at the barrier, owing to a larger deformation. 
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Fig. 3. Fission barrier in 238 U without (squares) and with (circles) the E c . m . (rel, shape) correction. 
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Fig. 4. The barrier (Q =125 b, left) and beyond the barrier (Q =200 b, right) configurations in 
238 U; density contours 0.01 fm -3 apart. 



4. Conclusions 

Correctly calculated energies of nuclear shapes with constriction become lower than 
in the standard approach, so such shapes, in particular scission configurations, are 
less excited with respect to more compact configurations. Here are some conse- 
quences of this correction for theoretical predictions: 

(i) Fission barriers with configurations close to scission are lowered by aproxi- 
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mately £ c . m . which weakly depends on the mass number A. This brings the 
calculated fission barriers in relatively light (A = 70 — 100) nuclei much closer 
to the data^. 

(ii) A smaller correction is expected for fission barriers in A w 200 nuclei, e.g. 
in 198 Hg (Fig. 1). This is consistent with very elongated, constricted shapes 
at the fission saddles in these nuclei, quite different, however, from scission 
configurations. 

(iii) Fission barriers with energy at the scission point close (within a few MeV) to 
that at the ground state become shorter, and this leads to a moderate decrease 
in fission half- lives (Fig. 3 for 238 U). 

Modifications are expected for barriers and half-lives for multipartition and 
multifragmentation; a scission configuration for triparition will be lowered by ~ 
2-Ec.m., that for the decay into four fragments by ~ 3-E c . m ., etc. One can also notice 
that the correction to b arri ers would tend to destabilize hypothethic hyperdeformed 
minima studied in Ref. 

Due to the magnitude of the correction, it lowers substantially fission barriers 
(up to the actinides) and modifies fission life-times, except for the superheavy nu- 
clei. Even there, the correction should be accounted for when considering fission 
dynamics. A related correction may be necessary in methods other than HF, unless 
they correctly and smoothly describe binding during fission and fusion. 
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